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ABSTRACT 

We present a study of the stellar and circumstellar properties of Class I sources 
using low-resolution (R~1000) near-infrared K- and L-band spectroscopy. We 
measure prominent spectral hues and features in 8 objects and use fits to standard 
star spectra to determine spectral types, visual extinctions, ii'-band excesses, 
and water ice optical depths. Four of the seven systems studied are close binary 
pairs; only one of these systems, Haro 6-10, was angularly resolvable. For certain 
stars some properties found in our analysis differ substantially from published 
values; we analyze the origin of these differences. We determine extinction to 
each source using three different methods and compare and discuss the resulting 
values. One hypothesis that we were testing, that extinction dominates over the 
ii'-band excess in obscuration of the stellar photospheric absorption lines, appears 
not to be true. Accretion luminosities and mass accretion rates calculated for 
our targets are highly uncertain, in part the result of our inexact knowledge of 
extinction. For the six targets we were able to place on an H-R diagram, our age 
estimates, <2 Myr, are somewhat younger than those from comparable studies. 
Our results underscore the value of low-resolution spectroscopy in the study of 
protostars and their environments; however, the optimal approach to the study 
of Class I sources likely involves a combination of high- and low-resolution near- 
infrared, mid-infrared, and millimeter wavelength observations. Accurate and 
precise measurements of extinction in Class I protostars will be key to improving 
our understanding of these objects. 
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Introduction 



Class I sources represent one of the earliest stages of star formation and are identi- 
fied by a rising spectral energy distribution (SED) at wavelengths longer than 2 /xm (Lada 
1987). These sources are deeply embedded within molecular clouds and are very faint or 
undetectable at visible wavelengths because of a thick envelope of circumstellar dust. This 
material effectively envelops the whole star and as a result, the light from these young stellar 
objects (YSOs) is absorbed and re-radiated in the infrared (IR). Therefore, most studies of 
Class I objects employ observation s in the near-IR or longer w a velength regimes where these 
sourc es are relatively bright (e.g., iGreene fc Ladalll996l . l2002l : iDoppmann et al.ll2005l : [Beck 
2007h . 



Following the suggestion of Lada (1987), it has been commonly thought that Class I ob- 
jects inhabit an earlier evolutionary stage relative to Class II sources (classical T Tauri stars, 
or CTTSs). However, showing that the central YSOs in Class I sources display signatures 
indicative of an earli er evolutionary phase compared to CTTSs has not been straightforward. 
Kenyon et al.l (119981 ) conducted a visible light, low spectral resolution survey of Class I YSOs 
in the Taurus star forming region (SFR). They identified preliminary spectral types and lu- 
minosities and determined that for Class I sources these properties were similar to those of 
the Class II objects in Taurus. However, these investigators did find a greater frequency 
and intensity of outflows a s deduced from forbidden hne emission in the Class I objects. 
White fc Hillenbrandl (120041 ) used high spectral resolution data in the visible to measure the 
stellar properties of several dozen Class I objects and reached similar conclusions, although 
they argue that the larger equivalent width values observed in forbidden emission lines from 
Class I sources might be attributable to the effect of circumste llar disk orientation causing 
obscuration of the cen tral continuum source ( White et al.ll2007l ) rather than to a higher in- 
cidence of strong jets. White fc Hillenbrand ( 20041 ) also found that the veihng and derived 
accretion r ates were similar for Cla ss I and Class II sources; however, using high-resolution 
IR spectra, iDoppmann et al.l (120051 ) found that the veiling and associated accretion rates of 
Class I and flat-spectrum YSOs are higher than those of the Class II objects. As a result 
of these discrepancies, it is unclear how much of the flnal mass of a YSO is accreted during 
the Class I phase. Episodic events have been pro posed to account for signiflcant growth over 
short periods (e.g., iKenyon et al.lll998l ): however, IWhite et al.l (120071 ) argue that the process 
by which protostars acquire the majority of their mass is still unconflrmed. 

From an evolutionary point of view, another property that is expected to differ between 
the two classes of YS Os is their rotation rat e s, but again unambiguous observational evidence 
for this is lacking. IWhite fc Hillenbrandl (120041 ) determined that their sample of Class I 
objects in Taurus is rotating at vsini <35 km s~^ on average, which makes them generally 
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indistinguishable from Class II sources in this SFR. However, ICovey et al.l (120051 ) found 
that Class I arid fiat -spectrum objects (apparently transitioning from Class I to Class II; 
Greene fc Ladalll997l ) in the Taurus and Ophiuchus SFRs do rotate more quickly than Class 
II objects in the same regions. To a ccount for rotational slowing between the Class I a nd Class 
II phases. iMontmerle et al.l (120001 ) invoked the onset of magnetic disk braking (e.g.. iKoenigl 
199ll : IShu et al.lll994l ) over the protostellar lifetime. In the framework of this paradigm, once 
disk locking has fully engaged, the star has evolved t o the Class II stage and its rotation has 
slowed accordingly. Alternatively, iMatt fc Pudrita (120051 ) proposed accretion driven winds 
to account for angular momentum changes in young stars. 

It seems unlikely that all Class I objects could be Class II systems with irregular appear- 
ances attributable to geometric effects. Some Class I sources no doubt are simply CTTSs 
seen through edge on disks, but the rest may well be part of an overall evolutionary sequence 
from Class O/I to Class I/II YSOs. Untangling the evolutionary state and geometric effects 
influencing these classifications is, however, apparently quite complex. 



Following the early, low-resolution (~600) near-IR spectroscopy studies of lGreene fc Lada 



(119961 . 1997) on Class I sources there ha s been a shift to primarily high spectral resolu 



tion ( R>18,000) studies of prqtostars (e.g 



2004 : jPoppmann et al.l l2005l : I Covey et al. 



2005 



Greene fc Ladall2000l . l2002l : IWhite fc Hillenbrand 
20061 ) as the requisite instrumentation has 



become available on large telescopes. High-resolution spectra permit the measurement of 
wsim, magnetic fields, and other detailed properties. Although this work has proved valu- 
able, the wide wavelength coverage available with low dispersion spectrographs provides 
unique leverage in the study of these very embedded young stars, allowing for the simul- 
taneous determination of the extinction and the examination and characterization of many 
spectral features. For objects in common with high dispersion studies, the results we obtain 
here, based on low-resolution IR data, are not always the same, suggesting that the two 
approaches are complementary. 

While a number of continuum and line emission processes likely play a role in shaping 
the emergent spectra of Class I objects, we undertook the current study to test the hypothesis 
that extinction is the dominant factor, as indicated by analysis of the protostar YLW 15A 
(jPrato et al.ll2003l ). In this current paper, we determine the properties of a small sample of 
8 Class I YSOs using a number of techniques such as measuring the equivalent width of key 
spectral absorption and emission line features, placing sources in a color-color diagram using 
photometric data obtained from the literature, and fitting our observed spectra to standard 
star spectra to determine effective temperature (Tg//), extinction, and K-hand excess. Our 
goal is to understand better the properties of protostars and the characteristics of their 
surrounding environments using photometry in combination with low-resolution spectroscopy 
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of atomi c and molec ular lines, as well as of solid state features such as water ice absorption at 
3.1 fim (iBeckl 120071 ). We examine several approaches to measuring extinction in our targets 
and compare these with each other and with estimates found in the literature. In this study 
we find a wide range in the strength and shape of protostar spectral features, both emission 
and absorption lines, from source to source, even in our relatively small sample of 8 objects. 
The principal unifying feature to the entire set is a rising spectrum across the ii'-band. The 
observations and data reduction procedures are described in §2. In §3 we present the analysis 
of the data. A discussion appears in §4 and the results are summarized in §5. 



2. Observations and Data Reduction 

Our sample consists of 8 Class I YSOs in the Taurus star-forming region, two of which 
comprise a relatively wide binary, Haro 6-10. Three of the YSOs were unresolved binaries, 
IRAS 04239+2456, L1551 IRS 5, and IRAS 04361+2547. The objects are listed in Table 1, 
along with coordinates and observing dates. The coordinates reported here were determined 
using the Two Micron All Sky Survey (2MASS) interactive image service. Initially, using 
coordinates in SIMBAD, 2MASS images were found to show a displacement of sometimes 10" 
or more between the coordinate center and the stellar image. By iterating the coordinates 
until these matched, we determined the best values for our sample target s. Haro 6-10 is 



unreso lved in 2MASS images. The sample was drawn from the study of I Greene fc Lada 



( ll996[ ) and selected on the basis of steeply rising K-hand spectra; Greene & Lada identified 
7 of the 8 systems as Class Is and IRAS 04489+3042 as a fiat spectrum source. In addition, 13 
spectral type standard stars were observed and are listed in Table 2 along with their spectral 
types, i^-band magnitudes, and dates of observation. These standards were selected to cover 
a wide range of low-mass spectral types for use in the model fitting described below. 

The data were obtained at the NASA Infrared Telescope Facility on Mauna Kea in 
2000 November 18— 20 (UT), using SpeX, the facility low- resolution near-IR spectrograph 



( iRayner et al.ll2003l ). A 1024x1024 InSb array was used to record the spectra. Simultaneous 
observations of the K- and L-bands (~2. 0—4.2 fim) were acquired by means of a prism cross- 
disperser. A slit width of 0'.'8 was used for all observations, yielding a spectral resolution 
of R = A / (5A = 1000. The seeing ranged from 0'.'5 to 0'.'9. The data were acquired in 
pairs, nodded along the slit between frames. Exposure times varied between 0.51 s and 30 
s, and multiple exposures were taken in each nod position. Stars of spectral type AO were 
also observed for removal of telluric features. Exposures of an argon lamp were taken for 
wavelength calibration. Flat field and dark exposures were also acquired. 
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Data were reduced using the REDSPEC codqj. Normalized flat fields were created by 
taking the median of the set of three flat field exposures taken for each star. The effects 
of subtracting median dark frames from the flat were insignificant so this differencing was 
not performed. Target exposures were medianed to create a single image of the spectrum of 
each star at each nod position. For each object these medianed nod pairs were differenced 
and divided by the corresponding flat. The images were spatially rectified using fourth- 
order fits to bright AO star traces in each order. When the REDSPEC code could not fit 
the traces automatically, as occurred in an L-band order where there was significant atmo- 
spheric contamination, manual fits to the traces were performed. Wavelength calibration 
was accomplished by fitting a second-order polynomial to the positions of identified lines 
in the argon lamp line spectrum. Several rows containing the stellar spectral data in the 
rectified, differenced image, one positive trace and one negative trace, were then summed. 
By subtracting the resulting negative spectrum from the positive, we accomplished a double- 
difference procedure that eliminated OH night sky emission line residuals. AO calibrator star 
spectra, observed at similar airmasses as the standard and target stars, were used for the 
removal of telluric absorption lines. Intrinsic atomic hydrogen lines were interpolated over 
in the calibrator star spectrum, which was then divided into the standard and target star 
spectra. The ratio was multiplied by a featureless 9500 K blackbody curve to restore the 
original shape of the continuum of the standard or target star. After reduction, cosmetic 
improvements (such as interpolation over bad pixels) were made to the spectra. The reduced 
/T-band spectra of the standard stars appear in Figures 1 and 2. The reduced K- and L-band 
spectra of the Class I objects appear in Figures 3 and 4. 

Four of our target systems are binaries; however, only Haro 6-10 was angularly re- 
solvable. When extracting the Haro 6-10 S and N spectra, the number of extracted rows 
was chosen so as to avoid contamination by the companion star. The spectra of IRAS 
04239+2436, L1551 IRS 5, and IRAS 04361+2547 are composed of the blended light of their 
primary and secondary components. 



3. Analysis 

3.1. Target Colors and Corresponding Extinction 

Near-IR magnitude s for the Class I samp le (Table 3) were obtained from the 2MASS 
point source catalogue (ISkrutskie et al.l l2006l ) and used to construct a J— H versus H— K 
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color-color diagram (Figure 5). Haro 6-10 N lacks a J-band measurement and hence is not 
shown. The location of Haro 6-10 S on the diagram implies strong contamination from 
circumstellar scattered light. All other targets show evidence for a large near-IR excess and 
most are highly extinguished. 

To measure the extinction from the color-color diagram, we used the relation = 
13.83(J - H)obs - 8.29{H - K)obs - 7.43, derived in IPrato et all J2003h . This dereddens 
the targets to the classica l T Tauri locus, i .e. the intrinsic colors of a young star plus 
circumstellar disk system (IMeyer et al.l 119971 ) . although it does not take into account the 
effects of disk geometry (e.g., flaring), disk inclination, or the accretion rate of disk material 
onto the star(s). Also, the classical T Tauri locus was determined for K and M stars only 
and thus is only strictly applicable to these spectral types. This process assumes that the 
envelope surrounding the Class I system contributes only reddening and that the intrinsic 
colors of the star-|-disk system are those of a classical T Tauri star. The extinction results 
appear in the last column of Table 3 and do not include Haro 6-10 N because of the missing 
J band data. 



3.2. Equivalent Widths and Surface Gravity 

Equivalent widths for spectral features in the i^'-band were measured for the standard 
stars and for the Class I YSOs and are given in Tables 4 and 5, respectively. The equivalent 
width of the Bra emission line for the Class 1 sources is also given in Table 5. Equivalent 
width uncertainties were determined by varying the location at which the relative continuum 
was measured. For the i^-band, this yielded uncertainties of ~0.5 A and for the L-band, ~2 
A. When no value is given in Tables 4 or 5, the equivalent width is less than the uncertainty 
or the line is not detected at all. 

Only two targets in our sample, L1551 IRS 5 and IRAS 04489-^3042, showed Nal, Cal, 
and CO (2—0) all in absorption. We c ompared the locations o f these two stars with the 



dwarf (from our data) and giant (from [Wallace fc Hinkld 119971 ) star Nal-|-Cal versus CO 



(2—0) equivalent width loci (similar to the loci data in Prato et al. 2003) and find that the 
Class I sources are consistent with giant star surface gravities. We conclude that, for Class 
1 protostars, at least in some cases, low surface gravity standards would be preferable as 
spectroscopic templates. 
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3.3. Spectral Types, K-hand Excess, and Extinction Revisited 

To determine the underlying spectral type, i^-band excess, tk, and extinction down 
to the protostellar photosphere, we compared the A'-band YSO spectra to a suite of dwarf, 
subgiant, and giant spe ctral type standard s (Figures 1 and 2). We initially followed the same 
procedure described in iPrato et al.l (120031 ) which we briefly review. For a range of trial Ay 



and tk values, the Class 1 target spectra are dereddened using the reddening law of Rieke 
and Lebofsky (1985) and the veiling continuum is removed. This modified spectrum is then 
fit to spectral type standard templates by determining the multiplicative scale factor that 
minimizes x^- The combination that gives the lowest overall value of then gives us an 
estimate of Ay, r^, and the spectral type approptiate for each Class I source. This procedure 
provided initial estimates of these key properties; however, additional experimentation with 
wavelength dependent veiling and visual examination of the resulting fits was used to arrive 
at the final values shown in Table 6. Crude uncertainties on the order of 5—15% were 
estimated from visual inspection of the YSO fits to the comparison standard stars. Some 
comments on the fitting process for each target follow. 



3.3.1. IRAS 04016+2610 



White fc Hillenbrand (120041) and lDoppmann et al.l (120051 ) identify a K star spectral type 



and lDoppmann et al.l (120051 ) find a large K-hernd excess for this protostar. Given the similar 
depths of the Cal, Mgl, and first CO bandhead evident in the spectrum of IRAS 04016+2610, 
an early K spectral type appears to provide the best match. The i^'-band excess has a positive 
slope of Fx/X ~0.5 and a y-intercept of 0.1. Unfortunately, the subtraction of such a large 
i^-band excess from the observed spectrum results in a noisier spectrum that the original; 
higher signal to noise data would improve the accuracy of the fit. Experimenting with 
multiple spectral type standard spectra, a range of extinctions from 43—53 mag, and veiling 
at 2.2 fim of 0.5—1.5, we conclude that IRAS 04016+2610 is best fit with a K3 spectral type. 



A„ = 45 mag, and 



1.3. This spectral type and excess are consistent with the results of 



Doppmann et al.l (120051 ) using high-resolution (R=18,000) i^T-band spectra. 



3.3.2. IRAS 04181+2654A 



The data of iBeckl (120071 ) taken of this source ~3 years after our spectra, and at higher 
signal to noise, show better-defined stellar absorption features. However, both the shape of 
the continuum and the depth of absorption lines, shown in Figure 3, provide important limits 
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from which to determine the best fit. The presence of water vapor in the at mospheres of M 



type objects produces a break in the slope of the spectrum at about 2.3 /im (IWilking et al. 



19991 ). Thus, objects with such a break (Figures 1 and 3) are unambiguously of later spectral 



type. In lieu of clearly detectable stellar absorption lines, we estimated an upper bound for 
the extinction by assuming a zero K-hand excess and fitting the slope of the continuum to 
an M star. The continuum spectrum of IRAS 04181+2654A is consistent with an early M 
star and an upper limit of 34 mag. iBeckl ( 20071 ) found an of 18 mag and an M3l'^^ 
spectral type. 



3.3.3. IRAS 04239+2436 

All spectral features, including Na I and CO, observed in the K band for this target 
are in emission. No break in the continuum is detected, thus a spectral type later than K7 
(Figure 1) is unlikely. Using a KO V standard and setting the ii'-band excess to zero we 
estimate an upper limit for the visual extinction of ~45 mag. 



3.3.4. Haro 6-10 S 



We found a constant i^-band excess of ~0.5 for Haro 6-10 S and an extinction of 
30 mag. No break in the continuum was observed, implying a spectral type earlier than 
MO. Small absorption lines of Nal, Cal, and Mgl are clearly visible; however, the CO is 
probably in emission, f illing in what would otherwise be fairly deep absorption features (see 



Doppmann et al.ll2008l ). We assign Haro 6-10 S an early K spectral type and use a K3 type 



to estimate an approximate temperature, although this is highly uncertain. 



3.3.5. Haro 6-10 N 

Setting the i^-band excess to zero we estimated an upper limit for the extinction and 
the spectral type of the emission line object Haro 6-10 N. We find an upper limit of 59 mag 
for the extinction and a spectral type consistent with an early K or G type star. 
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3.3.6. L1551 IRS 5 

Figure 6 shows our fit of L1551 IRS 5 to an M3 giant standard after acccounting for an 
extinction of 28 magnitudes and a constant excess of 1.0 across the i^'-band. 



3.3.7. IRAS 04361-2547 

With a i^-band excess of zero we find an upper hmit of 30 magnitudes for the extinction 
to IRAS 04361—2547; the curvature in the spectral continuum implies a spectral type later 
than M3. 



3.3.8. IRAS 04489+3042 

A hybrid standard star spectrum obtained by averaging the M3III and M4V spectral 
type standards provided the best fit to this source after accounting for a constant i^-band 
excess of 1.4 and 18 magnitudes of extinction. The deep CO bandheads and characteristic 
M star break at ~2.28 /xm support this conclusion. We find a negative slope of 0.5 for the 
excess across the K band and a y-intercept of 2.5. 



3.4. Ice Band Absorption and a Third Extinction Estimate 



The most obvious feature in the L-band spectra of Class I YSOs is the wide, deep, water 
ice absorption at 3.1 /im (Figure 4). This feature appears in all 8 of the target sources to 
varying degrees. In order to use the 3 /xm ice feature to estimate the optical depth through 
the envelope, it is necessary to isolate the ice absorption profile from the overall shape of 
the spectral energy distribution (SED) in the 2-4 /xm range. To do this, we construct a 
simple model which is fit to each source. The model consists of a template standard star to 
which we add continuuum veiling; this sum is then subjected to reddening. The reddened 
spectrum is then scaled by a multiplicative factor to account for the distance to and size of 
each source. The template star spectrum is our observed spectrum of either HR 753 (K3V) 
or GJ 876 (M4V), depending on which is closer in spectral type to o u r esti mate for each 
Class I source. The reddening law is again that of iRieke fc LebofskyI (119851 ). The veiling 
emission is taken to be quadr atic in wavelength. We then use the technique of Marquardt 
( iBevington &: RobinsonI Il992l ) to determine the best fitting parameters which match our 
model spectrum to the spectrum of each Class I source. In the fitting procedure, we fit only 
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the spectral regions shortward of 2.5 /im and between 3.65 and 4.1 /im. For Haro 6-10 N, 
we also ignore two strong absorption features, likely artifacts of the reduction, which appear 
at 3.76 and 3.88 /im. Once we have fit each observed spectrum, we divide that spectrum by 
the fit and set the minimum value of the ratio equal to e~'^^'"' and solve for Xjce- An example 
of our fit to IRAS 04016+2610 appears in Figure 7. 



Whittet et al.l (119881 ) found a linear correlati on between Xjce and the visual extinction 



Aj,, for the interstellar medium. ISato et al.l (119901 ) examined this correlation for a sample of 
protostars, assuming the ice to be in the circumstell ar material surrounding the star itself, 
and used the same relation as IWhittet et al.l (Il988l ). Xjce = 0.093(74y — Avc), adopting a 
critical visual extinction of Aye = 1.6 mag. We use this modified relation to calculate the 
extinction as traced by the circumstellar ice feature. The results appear in Table 7. 

The combination of near-simultaneous moderate resolution spectroscopy from the near 
to the mid-IR allows for the comparison of solid state features and hence reveals physical and 
compositional properties of Class I envelopes (e.g., iFurlan et al.l 120081). We examined the 



correl ation of the 10 fim silicate absorption feature strength, Fx^ from iKessler-Silacci et al. 



(120051 ). calculated from a Gaussian fit to the normalized mid-IR spectra, with optical depths 
we determined for the 3.1 fim water ice feature (Table 7). Figure 8 shows a trend towards 
higher ice optical depths for strong silicate absorption sources only; no significant correlation 
is obvious. 



3.5. H-R Diagram 



To compare the relative masses and ages of our Class I sources with those of Class II 
objects, we plotted stars for which we have the most reliable spectral types, determined 
from both absorption lines and continuum sh ape, on the H-R di agram (Figure 9). For 



four targets we used luminosity estimates from iFur 
04016+2610 and Haro 6-10 S, we used values from 



an et al. (2008). For two others, IRAS 



Doppmann et al.l (120051 ) either because 



the value from iFurlan et al.l (12008) was anoma l ously low (implying a main sequence age 



for IRAS 04016+2610) or because IFurlan et al.l (120081 ) did not observe the source (Haro 6- 
10 S). We used o ur spectral type estimates (Table 6) combined with the dwarf conversion 
(Table 2) given in I JohnsonI (Il966l ). to determine Teff- Interestingly, J ohnson's spect r al typ e 



to Teff conversions closely parallel those derived for young stars by iLuhman et al.l (120031 ). 
from G through mid-M spectral typ es. In Table 8 we li st our Tpj-y, and f or co mparison 
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were possible to derive. 

Of the seven systems studied in this paper (taking Haro 6-10 S and N as one system), four 
are close binaries (Table 1). The three single systems in our sample along with Haro 6-10 S 
and two of the unresolved binaries appear in Figure 9. There is an inherent overestimate of 
the stellar luminosity for the unresolved pairs, which are treated as single objects by Furlan et 
al. (2007). Also, geometrical effects produced by the relative orientations of the circumstellar 
disks of each binary component could lead to anomalous measurements of extinction and 
excesses, for example, in the case of a configuration in which one star is obscured by the 
disk of the other. To the extent that circumstellar disks (flared or seen nearly edge on) 
contribute to the extinction, it is quite possible the extinction is very different to the two 
members of the binary. Without angularly resolved observations of the component objects in 
these systems, the degree to which binarity distorts the results of Class I studies, specifically 
the effective temperature and luminosity estimates and thus a source's location in the H-R 
diagram, is unknown. The unresolved binaries are located far above the evolutionary model 
tracks in Figure 9, suggesting unrealistically young ages for these systems. 



3.6. Accretion Luminosity and Mass Accretion Rate 



For comparison with other studies of Class I objects as well as with Class II targets, we 
calculated the accretion luminosity. Lace, of the eight stars in our sample using the correlation 
between Lace and Br7 line luminosity found by Muzerolle et al. (1998). Our calculation is 
based on the i^-band magnitude, the Br7 emission line strength, and the overall extinction, 
for which we used the value determined by model fitting to standard stars as described in 
§3.3. The results appear in column 2 of Table 9. In most cases, the measurement is significant 
only at the one sigma level because the inherent un certainties in the relat ionship between the 
accretion luminosity and the Br7 line luminosity (IMuzeroUe et al.lll998l ) are large. For half 
our sample we use upper limits for the extinction; if instead we use the values for extinction 
determined from the ice feature (Table 7) we obtain lower values for Laec in the cases of 
IRAS 04181+2654A {log{Laec)= -0.84±1.12), IRAS 04239+2436 (\og{Laec)= -0.19±1.06), 
and Haro 6-10 N {\og{Laec)= -0.94±1.14). For IRAS 04361-2547, the other system with 
an upper limit only for Ay in Table 6, the ice feature (24.2 mag) is similar to the upper 
limit (<30 mag). 

Fo r the six targets in our sample for which we have estimates of both Tg// and the lumi- 
nosity (IFurlan et al.ll2008f). we calc u lated the stellar radius, R*, and hence the mass accretion 
rate, Ivl. following ICuUbring et~aD Jl998h : M = 1.25LaecR*/{GM^) (Table 9). Stellar mass 
(M*) was estimated from the location of the targets in the H-R diagram. Primarily as a re- 
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suit of the very large uncertainties in the calculated accretion luminosity, Lace, uncertainties 
in our mass accretion rates are all approximately an order of magnitude. 



4. Discussion 
4.1. Comparison with Previous Results from Literature 

4-1.1. Extinction 



As demonstrated here and also shown in iBeckl (120071 ) , different approaches to measuring 



the extinction along the line of sight to Class I protostellar photospheres produce widely 
different results. We have used the position of sources in a near-IR color-color diagram, 
fitting to absorption features and the continuum shape of the i^'-band spectra, and fitting to 
the L-band water ice feature to determine and have found values that differ significantly. 
Figure 10 shows the extinctions determined from de-reddening objects plotted on the color- 
color diagram and from i^-band spectral fitting as a function of determined from fitting 
the L-band water ice feature (Figure 7). Color-color diagram determined extinctions appear 
to scatter randomly, whereas extinctions measured from our i^-band spectra, excluding the 
upper limits, show some positive correlation with the ice A^. This suggests that the color- 
color diagram approach does not trace the same material along the line of sight. The color- 
color approach was based on the possibly flawed premise that the circumstellar envelope 
present in a Class I system contributes only reddening and assumes that the intrinsic colors 
of Class I protostars are the same as those of classical T Tauri stars. As a counterexample, 
consider objects that fall below the classical T Tauri star locus in the color-color diagram, 
i.e. IRAS 04016+2610 and Haro 6-10 S. These are no doubt extinguished; however, the 
effect of scattered light in the surrounding circumstellar environment yields unusual colors. 
Measurements of the extinction from the optical depth of the water ice feature in the L-band 
probes the extinction through the region where this feature forms, while possibly missing 
extinction caused by dust in inner, warmer regions of an envelope or disk where it is too hot 
for water ice to survive. This could explain the systematically higher A^ values from the 
i^'-band spectroscopic fit approach (again, excluding the upper limits) which show a similar 
slope as the one-to-one correlation and presumably take into account all sources of extinction 
along the entire line of sight to the target system. 



C ompa rison of our extinction results with those of IWhite fc Hillenbrand! (120041 ) and 



Beckl (120071 ) reveals that these studies are generally consistent with our smaller estimates 



of ext inction based on fitting of the L-band ice absorption feature. IWhite fc Hillenbrand 



( I2OO4J ) provide extinctions determined from J— H colors and iBeckl (120071 ) derives extinctions 
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using several different approaclies, inclu ding a spectral fitt i ng pr ocedure similar to the K- 
band fits employed here. In their §3.8, iDoppmann et al.l (120051 ) describe an approach to 
estimating extinction, based on /T-band magnitudes, in order to derive stellar luminosities; 
however, they do not provide values of extinction in their paper. It is difficult to know which 
approach yields the "correct" extinction to a protostar. Ideally, a combination of near- 
simultaneous high- and low-resolution spectroscopy could provide careful measurements of 
underlying absorption line ratios as well as a measurement of the continuum slope. Even 
relatively small uncertainties in extinction can have a strong impact on derived Class I object 
properties, particularly accretion rates, and hence studies of these targets can be significantly 
improved with more attention to this problem. 



4- 1-2. Spectral Types, Surface Gravity, and K -band Excesses 



Three of our targets, IRAS 04016+2610, Haro 6-10 S, and IRAS 04489+3042, fo r 
which we deterrn i ned T e ff^ and r^, w ere also analvzed bvlWhite fc Hillenbrand! (120041 ). 
Doppmann et al. ( 2005h. Luhman ( 2006), a nd Doppmann et al. ( 20081 ). Another, L1551 IRS 
5, was analyzed by lDoppmann et al.l (120051 ). Our results for spectral types, surface gravities, 
and /T-band excesses typically agree to within l—2a with those from the literature with the 
except ion of Haro 6-10 S and L1551 IRS 5, discussed below. For IRAS 04489+3042. iLuhman 
(120061 ) found a spectral type of M3.5— M4.5 on the basis of visible light observations, and 
M3— M4 on the basis of IR observations, in excellent agreement with our determination of 
an M3.5 spectral type for this source. 



For Haro 6-10 S, iDoppmann et al.l (120081 ) used high-resolution, high-signal to noise K- 
band spectra to derive a Tg// of 3800 K, cooler than our value of ~4800 K, a surface gravity 
(log g) of 4.0, and a iT-band excess of 2.5, substantially larger than our uncertain estimate of 
0.5. Based on the appearance of their spectra, these values for Haro 6-10 S are more reliable 
than our findings and result in a mass estimate smaller than ours by a factor of >3; both 
studies result in similar ages (2—3 Myr). 



For L1551 IRS 5, we find an excellent match to an M3 HI spectral type. IDoppmann et al. 



(120051 ) find a Te// of 4800 K, corresponding to an early K spectral type, and a log g of 4.0. 
In both our analysis and theirs, the i^-band veiling is ~1. Figure 6 shows our M3 HI 
fit of the modified young star spectrum, with =28 m ag and a constant i^-ba nd excess 
taken into account. Inspection of the spectra presented in IDoppmann et al.l (120051 ) shows an 
inconsistency with their earlier spectral type and higher gravity, possibly as the result of a 
degeneracy between the equivalent widths of the Na I and Mg I/Al I lines as a function of 
Tg// and log g (Doppmann 2007, priv. comm.). The lower surface gravity found in our fit 
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is also consistent with our comparison of the location of L1551 IRS 5 on a plot of the Na I 
+ Ca I versus CO (2-0) equivalent widths (§3.2): it clearly lies along the giant star locus. 

In general, from our i^'-band spectral fits we find higher extinctions and lower excesses 
than those derived in previous Class I studies. We have used arbitrary slopes for fitting 
the i^-band excess, with no physical basis in a realistic model of circumstellar disk and/or 
envelope radiation. A combination of high signal-to-noise ratio low- and high- resolutio n 
spectra and detailed modeling of the spectral energy distribution, as in lFurlan et al.l (120081 ) . 
will be necessary to improve these uncertain results. 



4.1.3. Accretion Luminosities and Mass Accretion Rates 



For six of our tar gets, published ac cretion luminosities are available for comparison in 
MuzeroUe et al.l (119981 ) and iBeckl (120071 ): these are provided in Table 9. Although the pro- 
cedure for measuring the r aass accretio n rates followed in this paper and in the work of 
MuzeroUe et al.l (119981 ) and iBeckl (120071 ) was the same, i.e. based on the Br7 emission line, 
our values are all larger by >1— 2 a. Four objects have upper limits only for the spectro- 
scopically determined A^. We recalculated the accretion luminosities for the three largest 
upper limits using the A^ values based on the 3.1 /im ic e feature (see ^3 . 6). T he recalcu- 
lated Lace values are much closer to the values found by iMuzeroUe et al.l (119981 ) and iBeck 



(120071 ) for the same targets, illustrating that the relatively large values and upper limits for 
the extinction determined from our i^-band spectra are responsible for our relatively large 
accretion luminosities (Table 6). The impact of these larger accretion luminosities propa- 
gates into our calculation of the mass accretion rates, which range from ~ 5 x 10~^ Mq 
yr~^ to ~ 2 X 10^^ Mq yr~^ (Table 9), calculable for six of our targets. For the targets 
with published accretion rates, our estimates are equivalent or larger by up to two orders of 
magnitude (Table 9), however, with an uncertainty of approximately an order of magnitude. 
Given this large uncertainty, our small sample size, and the ambiguity in the determination 
of A^, it is impossible to draw definitive conclusions regarding potential enhanced accre- 
tion in Class I versus Class I I objects. Class II mass accretion rates are typically >10~^ 
(iHerczeg fc Hillenbrandll2008l ). Our results point to larger accretion luminosities and mass 
accretion rates than those of Class II objects, sugges ting an evolutionary pro gression from 
protostars to classical T Tauri stars, as discussed in iDoppmann et al.l (120051 ) for example, 
but without the signal to noise to support this conclusion at a significant level. 
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4.2. Accreting Protostars or Misclassified T Tauris? 



Recent st udies of Class I protos t ars d iffer in their conclusions regarding the nature of 
these objects. IWhite fc Hillenbrandl (120041 ) suggest that most Class Is have moved beyond 
their main accretion phase and might actua l ly be Class II stars seen edge on through cir- 



cumstellar accretion disks. iDoppmann et al.l (120051 ) conclude that Class I objects are indeed 



actively accreting protostars, albeit spanning a range of accretion activity. Since the publica- 



tion o f both papers, their authors have reached some consensus and conclude in IWhite et al. 



(120071 ) that, at least for the Taurus SFR, one-third to one-half of the Class I objec t s are likely 
to be misclassified Class II stars seen through optically thick disks. I White et al.l (120071 ) also 
conclude that most Class I sources posses disk accretion rates below the expected envelope 
infall rates. Misclassified T Tauri stars would certainly show low disk accretion rates; how- 
ever, suc h objects would typic ally be sub-luminous by a factor of five compared to bona fide 
Class Is dwhitnev etaPboOsh . 



Another possibility raised in IWhite et al.l (120071 ) is that our knowledge of disk accretion 
rates and mass infall rates is incomplete. The uncertainties in these calculations are signif- 
icant: in this paper as well as in the literature these values are only known to a precision 
of about an order of magnitude. Our finding of relatively large accretion rates (Table 9) 
compared to previously derived values for Class I objects is not significant. Objects that 
are found to be highly extinguished, such as our results suggest, naturally have large accre- 
tion luminosities. However, the large uncertainties associated with the mass accretion rate 
estimates mean that we cannot be confident that our results indicate that these protostars 
are currently experiencing a particularly active accretion phase. Our targets do, however, 
lie relatively far up on the evolutionary tracks in the H-R diagram (Figure 9), with ages 
<2 Myr, consistent with a population younger than that of the Class II T Tauri stars with 
ages of ~1 to several Myr. 

Another source of uncertainty in the determination of Class I protostar properties is 
their variability. Class I YSOs show consid erable variation in their spectra, fluxes, colors, 
and excess IR emission. iLeinert et al.l (120011 ) studied the variability of the Haro 6-10 system 
over a twelve year period and found substantial change in the magnitudes and colors of the 
two components, especially at shorter wavelengths, i^'-band veiling measurements for the 
Class I protostar YLW 1 5A (IRS 43) varied significantly over ^5 years, possibly on tina e 
scales as short as weeks (ILuhman fc Riekd Il999l : iGreene &: Ladal l2002l : jPrato et al.l |2003| ) . 
These results suggest that the geometry of the obscuring material in front of the protostellar 
surface may be shifting rapidl y, perhaps producing variations in the extinction along the line 
of sight to the protostar (e.g., iLeinert et al.ll200ll : iBeck et al.ll200ll ). 



Consistent with a younger evolutionary stage for the Class I sources is the ubiquitous 
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atomic and molecular hydrogen emission in all but one of our targets (Figure 3 and Dopp- 
mann et al. 2005). The atomic hydrogen emission is likely associated with disk accretion 
onto the central star and the molecular emission with shocks created by infalling or out- 
flowing gas interacting with the circumstellar environment. IRAS 04489+3042 shows only a 
small Br7 emission line and does not reveal any molecular hydrogen emission, either in our 
study or in that of iDoppmann et al.l (120051 ) . All other targets have the H2 v = 1 — S(l) 
line in emission and about half also show emission in the H2 f = 1 — S(0) line. 



5. Summary 



We obtained K- and L-band spectra of a sample of eight Class I sources in the Taurus 
SFR with the SpeX spectrometer on the NASA IRTF 3 m telescope. We measured absorption 
and emission line equivalent widths, determined the optical depth of the L-band water ice 
feature, and, where possible, estimated the spectral type, visual extinction, and K-hand 
excesses for each of our targets. Our analysis and results are summarized below. 

- Using our derived spectral types in combination with 2MASS photometry and lumi- 
nosity estimates available ir i the literature, A ve pla ced our targets in the H-R diagram and 
used evolutionary tracks of iPalla fc Stahlerl (119991 ) to estimate stellar masses and ages. In 
addition, we calculated accretion luminosities from the Br7 line luminosity and determined 
the associated mass accretion rate for our targets. The six objects which we were able to 
place on the H-R diagram show ages of <2 Myr and clump into two groups with masses 
~0.2 Mq and 1.9 Mq. 

- For two objects we obtain significantly different results from those presented in the 
literature, Haro 6-10 S and L1551 IRS 5. The majority of ou r targets are multiples; al- 
though we resolved Haro 6-10 S and N, iDoppmann et al.l (120081 ) suggest that the southern 
co mponent is a spe c trosco pic binary. The high-resolution, very high signal to noise data 
of IDoppmann et al. J2008h indicate a later spectral type, closer to MO, than the early K 
type found by us . The companion to L1551 IRS 5 is seen only at very long wavelengths 
(>submillimeter; iLooney et al.l 119971 ) and it is thus unlikely that this component contami- 
nates the primary spectrum in our analysis. However, we found a robust fit for this target 
to an M3III spectr al type standard, a signi ficantly lower gravity and later spectral type than 
the ~K3V type of IDoppmann et al.l (120051 ). 



- Our experimental hypothesis, that obscuration of the photospheric absorption lines 
in Class I targets is dominated by extinction local to the sources, is not correct. Veiling of 
lines from the near-IR excess also plays a key role. It is likely that emission, for example 
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from warm CO in a circumstellar disk, also fills in spectral absorption lines in some cases, 
hampering the determination of the underlying protostellar characteristics. 

- We compared our derived properties with values from the literature and find that our 
larger estimates of extinction lead to larger accretion luminosities and larger mass accre- 
tion rates. We stress that the very large uncertainties, inherent in the relati onship between 



emission line luminosities and accretion luminosities ((Muzerolle et al\\l99ai) as well as in 
our observations, render these estimates highly uncertain, not only in our study but also in 
numerous examples in the literature. 

- The higher mass accretion rates we find are more consistent with the interpretation 
that Class I protostars are undergoing a relatively active mass accretion phase; however, 
because of the large uncertainties in these results, additional observations are required to 
confirm this interpretation. High-resolution, high signal to noise IR spectroscopy provides 
a promising approach, however, low-resolution data across the near-IR may be particularly 
useful for the evaluation of extinction. 

- The visual extinction to protostars appears to be an extremely important and poorly 
determined property, impacting estimates of all the key stellar properties including lumi- 
nosity, mass, age, accretion luminosity, and mass accretion rate. It therefore represents a 
critical area for improvement. 
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Table 1. Sample 



Object 


a 


5 UT Date of 


Name 


(J2000.0) 


(J2000.0) Observation 



IRAS 04016+2610 


04 


04 


43.0 


+26 


18 


57 


2000 


Nov 


18 


IRAS 04181+2654A 


04 


21 


11.5 


+27 


01 


09 


2000 


Nov 


18 


IRAS 04239+2436'^ 


04 


26 


57.1 


+24 


43 


36 


2000 


Nov 


18 


Haro 6-10^ 


04 


29 


24.2 


+24 


33 


00 


2000 


Nov 


19 


L1551 IRS 5^ 


04 


31 


34.0 


+ 18 


08 


05 


2000 


Nov 


18 


IRAS 04361+2547^^ 


04 


39 


13.9 


+25 


53 


21 


2000 


Nov 


19 


IRAS 04489+3042 


04 


52 


06.7 


+30 


47 


18 


2000 


Nov 


20 



^Binary of separation 0'.'30 (jReipurth et al.ll2000l ). 



^Binary of separation 1'.'2 (ILeinert fc Haad 119891 ): listed as 
Haro 6-10 S (visible light primary) and N in subsequent tables. 
Haro 6-10 S is a rad ial velocity variable an d thus possibly a 
spectroscopic binary (jPoppmann et al.ll2008l ). 



^Binary of separation 0'.'35 (ILooney et al.l 119971 ) . 
'^Binary of separation 0'.'31 (ITerebey et al.lll998l ). 
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Table 2. Spectral Type Standards 



Object 


Spectral 


Ks 


UT Date of 


Name 


Type 


(mag) 


Observation 


HR 996 


G5 V 


3.0 


2000 


Nov 19 


HR 995 


G6 IV 


4.1 


2000 


Nov 20 


HR 7957 


KO IV 


1.4 


2000 


Nov 20 


HR 166 


KO V 


4.0 


2000 


Nov 19 


HR 753 


K3 V 


3.5 


2000 


Nov 19 


GL 846 


M0.5 V 


5.3 


2000 


Nov 18 


GL 908 


Ml V 


5.0 


2000 


Nov 18 


GL 15A 


Ml. 5 V 


4.0 


2000 


Nov 20 


GL 806 


M2 V 


6.5 


2000 


Nov 19 


GL 752A 


M3 V 


4.7 


2000 


Nov 18 


HR 9064 


M3 HI 


-0.1 


2000 


Nov 20 


GL 876 


M4 V 


5.0 


2000 


Nov 20 


GL 83.1 


M4.5 V 


6.6 


2000 


Nov 18 



Table 3. Sample Photometry 



Object 


J 


H 




L 


L-band 


Color- Color 


Name 


(mag) 


(mag) 


(mag) 


(mag) 


References 


Diagram A„ (mag) 


IRAS 04016+2610 


14.01+0.07 


12.16+0.08 


9.84^^ 


6.8^ 


1 


+1.6 


IRAS 04181+2654A 


16.22+0.08 


12.65+0.02 


10.34+0.03 


8.5*^ 


2 


22.8 +1.2 


IRAS 04239+2436 


15.75+0.09 


12.35+0.04 


9.99+0.02 


7.2'^ 


2 


20.0 +1.4 


Haro 6-10 S'^ 


11.54+0.03^ 


10.6+0.1 


8.6+O.P 


6.1+0.1^ 


3 


+1.9 


Haro 6-10 




13.5+0.1 


8.7+O.P 


4.9+0.1^ 


3 




L1551 IRS 5 


13.71+0.06 


11.51+0.05 


9.82+0.04 


7.4^ 


4 


9.1 + 1.2 


IRAS 04361+2547 


16.44+0.12 


13.02+0.04 


10.72+0.03 


8.9^^ 


5 


20.8 + 1.8 


IRAS 04489+3042 


14.43+0.03 


12.02+0.02 


10.38+0.02 


1.5^ 


6 


12.3 + 0.6 



Note. — Data are from 2MASS unless otherwise specified. 
^No uncertainty available. 
^L' magnitude. 

'^H, K, and L' magnitudes quoted for the highly variable Haro 6-10 system are from reference (3) and 
correspond to the date of their observations closest in time to our spectroscopy. 

•^AU J-band flux from Haro 6-10 is assumed to come from the "primary", southern component; the J 



magnitude is from 2MASS and was not observed concurrently with the H, K, and L data. 



References. — (1) Benson et al. (1984), (2) Beck (2007), (3) Leinert et al. (2001), (4) Cohen & Schwartz 
(1983), (5) Kenyon et al. (1990), (6) Myers et al. (1987) 



Table 4. Standard Star X-band Equivalent Widths (A) 





Spectral 


2.166//m 


2.208/im 


2.264//m 


2.281//m 


2.294/xm 


2.323/xm 


Name 


Type 


H I 


Nal 


Gal 


Mgl 


GO(2-0) 


G0(3-l) 


HR 996 


G5 V 


3.2 


1.3 


1.2 


1.1 


1.5 


2.1 


HR 995 


G6 IV 


2.2 


0.9 


1.1 


0.7 


4.7 


5.0 


HR 7957 


KO IV 


2.1 


0.9 


1.4 




5.9 


7.1 


HR 166 


KO V 


2.1 


1.5 


1.6 


1.7 


4.5 


6.5 


HR 753 


K3 V 


1.4 


2.2 


2.1 


1.9 


6.6 


6.6 


GL 846 


M0.5 V 




4.4 


4.3 


1.1 


7.4 


6.2 


GL 908 


Ml V 




1.9 


2.7 




4.8 


5.1 


GJ 15A 


M1.5 V 




3.3 


3.2 




5.6 


5.3 


GL 806 


M2 V 




3.8 


3.9 




6.6 


6.3 


GL 752A 


M3 V 




4.6 


4.5 




7.4 


7.1 


HR 9064 


M3 HI 




2.2 


2.8 


2.0 


20.9 


13.7 


GJ 876 


M4 V 




6.2 


4.2 




4.5 


6.2 


GL 83.1 


M4.5 V 




5.1 


2.4 




6.0 


4.3 



Table 5. Target Star Equivalent Widths (A) 





2.059/^m 


2.122^m 


2.166jMn 


2.208|[im 


2.264pm 


2.281jjm 


2.294pm 


2.323pm 


4.052pm 


Name 


He I 


H2 


H I 


Na I 


Ca I 


Mg I 


CO(2-0) 


CO(3-l) 


Br a 


IRAS 04016+2610 




-2.2 


-2.3 


1.0 


1.0 


2.0 


2.3 




-11.5 


IRAS 04181+2654A 


-1.4 


-1.9 


-9.9 


1.3 


1.6 


-0.8 




4.0 


-23.0 


IRAS 04239+2436 


-2.6 


-2.3 


-15.3 


-2.4 


1.4 




-11.9 


-10.5 


-26.0 


Haro 6-10 S 




-2.1 


-4.1 




1.0 


1.0 






-33.0 


Haro 6-10 N 




-2.5 


-3.1 


-1.1 






-6.7 


-2.2 




L1551 IRS 5 




-2.5 


-1.7 


2.4 


1.9 




16.7 


14.8 


-4.5 


IRAS 04361+2547 




-2.8 


-2.6 




0.6 


1.4 






-2.0 


IRAS 04489+3042 


-2.6 




-2.7 


2.0 


1.7 




4.6 


4.5 
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Table 6. Target Star Results 



Name 


Spectral 
Type 


Ay 








Log(L„cc) 
Lo 


IRAS 04016+2610 


K3 V ± 2^ 


45 ± 2 


1.3 


± 


0.2 


0.54 ± 0.99 


IRAS 04181+2654A 


Ml - M4 


<34 








0.35 ± 1.04 


IRAS 04239+2436 


G - early K? 


<45 








1.36 + 0.95 


Haro 6-10 S 


early K 


30 + 2 


0.5 


+ 


0.1 


0.37 + 1.04 


Haro 6-10 N 


G - early K? 


<59 








2.06 + 0.91 


L1551 IRS 5 


M3 III ±1 


28 + 2 


1.0 


+ 


0.2 


-0.69 + 1.12 


IRAS 04361+2547 


M4 V ± 2 


<30 








-0.80 + 1.13 


IRAS 04489+3042 


M3.5 IV ± 2 


18 ± 1 


1.4 


± 


0.4 


-1.26 ± 1.17 



^Uncertainties in spectral type refer to number of spectral subclasses. 
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Table 7. Ice Feature Optical Depths 



Name 






IRAS 04016+2610 


2.8 


31.7 


IRAS 04181+2654A 


1.0 


12.4 


IRAS 04239+2436 


1.4 


16.7 


Haro 6-10 S 


0.7 


9.1 


Haro 6-10 N 


0.5 


7.0 


L1551 IRS 5 


1.7 


19.9 


IRAS 04361+2547 


2.1 


24.2 


IRAS 04489+3042 


0.7 


9.1 



Note. — Uncertainties in Tjce are 
~0.1 and in are 1.1 mag. 



Table 8. Input and Results from H-R Diagram 







L. (L©) 


(Lq) 




L. (L©) 


Ago" 


Mass'' 


Name 


(K) 


Furlan et al. (2008) 


White & Hillenbrand 


f20041 


Doppmann et al. (2005) 


(years) 


(Mq) 


IRAS 04016+2610 


4800 ± 200 


0.45 ± 0.23 


0.45 




4.9 


I06'' 


2.0'^ 


IRAS 04181+2654A 


3350 ± 100 


0.14 ± 0.04 








2 


0.2 


IRAS 04239+2436 




0.15 ± 0.51 












Haro 6-10 S 


4800 ± 200 




1.8 




3.3 


206'' 


l.S" 


Haro 6-10 N 
















L1551 IRS 5 


3300 ± 200 


7.50 ± 1.46 






2.6 


<lo6 


0.2 


IRAS 04361+2547 


3200 ± 200 


3.60 ± 0.41 








<1g6 


0.2 


IRAS 04489+3042 


3250 ± 200 


0.21 ± 0.03 


0.2 




0.3 


I06 


0.2 



^Agc uncertainties are --^1 Myr. 
'^Mass uncertainties are 0.1 — 0.2 Mq . 
Based on luminosities from [p^oppmann et al ] i2005l ): all other results based on luminosities from [Furlan et al ] i2008l ). 



Table 9. Accretion Luminosity and Mass Accretion Rate 





log{L„ce/L0) 


log{Lacc/L0) 


log{Loce/L0) 


M (Mq) 


M (Mq) 


Name 


this work 


Muzerolle et al. (1998) 


Beck f 20071 


this work 


other 


IRAS 04016+2610 


0.54 ± 0.99 


-1.52 




2.20-7 


7.10-8" 


IRAS 04181+2654A 


0.35 ± 1.04 




-1.22 


4.90-7 


9.00-9*= 


IRAS 04239+2436 


1.36 ± 0.95 


-0.62 


-0.15 






Haro 6-10 S 


0.37 ± 1.04 






1.40-7 


2.Do-7'» 


Haro 6-10 N 


2.06 ± 0.91 


-1.34 








L1551 IRS 5 


-0.69 ± 1.12 






3.4e-7 




IRAS 04361+2547 


-0.80 ± 1.13 


-1.10 


-0.70 


1.90-7 




IRAS 04489+3042 


-1.26 ± 1.17 


-1.70 




1.6G-8 


4.10-10" 



IWhite fc Hillenbraiidl )2004l) . 
' iBeckl 120071) . 
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2.1 2.2 2.3 2.4 

Wavelength {/am) 

Fig. 1. — K-hand spectra of the dwarf and subgiant spectral type standards. The stars are 
identified in Table 2. Prominent spectral features are indicated. The spectra were normalized 
to unity at 2.2 /xm and shifted by a constant for presentation. 



-32- 




Fig. 2. — i^-band spectrum of the giant spectral type standard. The star is identified in 
Table 2. Prominent spectral features are indicated. The spectrum was normalized to unity 
at 2.2 /im and shifted by a constant for presentation. 
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T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

CO Ai/=2 

Hj H No Co Mg Bondheods 




2.1 2,2 2.3 

Wavelength (/xm) 

Fig. 3. — i^-band spectra for our sample of Class I sources. Prominent spectral features are 
indicated. The spectra were normalized to unity at 2.2 //m and shifted by a constant for 
presentation. 
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2.5 3.0 3.5 4.0 

Wavelength (/xm) 

Fig. 4. — K- and L-band spectra for the spectroscopic sample; the region around 2.7 //m is 
contaminated by terrestrial water absorption is not shown. The spectra were extracted from 
five overlapping orders and scaled to match in the overlap regions. The overall spectrum 
was normalized at 2.25 ^xm and shifted by a constant for presentation. The ice band feature 
and hydrogen lines are indicated. Pfund series emission lines are contaminated in part by 
residual inherent AO star absorption lines (Brackett lines are not contaminated). 
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2.5 3.0 3.5 4.0 

Wavelength (/xm) 

Fig. 4.— Con'd. 
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Fig. 4.— Con'd. 
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IRAS 04361+2547 / 



• IRAS 04181+2655 



IRAS 04239+2436 



L 1551 IRS 5 



• IRAS 04489+3042 



IRAS 04016+2610 



Haro 6-10 S 




t 

CTTs Locus 



I ^ I 




Ay^lO mag 



t—l I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I L. 



12 3 4 

H - K (mag) 

Fig. 5. — J-H, H-K color-color diagram of the Class Isample, except for Haro 6-10 N. Mag- 
nitude data were obtained from 2MASS and from the hterature. The dotted line separates 
objects with (to the right and below) and without a near-IR excess. The dwarf (dash-dot) 
and giant (sohd) star loci are overplotted at the bottom left. The CTTs locus and Ay — 10 
mag reddening vector from Meyer et al. (1997) are also shown. 
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2.1 2.2 2.5 

Wavelength (yC^m) 



2.4 



Fig. 6. — The ii'-band spectrum of L1551 IRS 5 (grey), accounting for 28 magnitudes of 
visual extinction and r^^ = 1.0, overplotted with the best fit stellar template (black), M3 III. 
The dotted line s h ows a K3V standard star, corresponding to the spectral type identified by 
Doppmann et al.l (120051 ). 
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2.0 2.5 5.0 5.5 4.0 

Wavelength {fim) 



Fig. 7.— The full K- and L-band spectrum of IRAS 04016+2610 (grey) and the best fit 
modeled stellar template (black) used to determine the optical depth of the 3.1 fim water 
ice feature. The broad, noisy section of spectrum from 2.6—2.8 /xm corresponds to a region 
of low transmission in the Earth's atmosphere. 
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2.5 - 



2.0 - 



^ 1.5 



1.0 



0.5 



0.0 



I I I I I I I I I I I I I 



-0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

Fx. 



Fig. 8. — The 3.1 /im ice optica l depth compared to the s trength of the 10 fim sihcate 
absorption feature as measured by iKessler-Silacci et ahl (120051 ) . The data show a weak trend 
towards higher values of Tjce among objects with the strongest silicate features. 
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Fig. 9. — Class I targets, for which we have estimated Tp f f and use values for the luminosity 
from iFurlan et al.l (120081 ) (circles) and iDoppmann et al.l (12005) fdia monds). plotted on an 
H-R diagram with the evolutionary tracks of iPalla fc Stahlerl (jl999l ). At the low- mass end, 
only every other mass track is labeled. Isochrones are plotted for ages of le6, 3e6, 5e6, 
le7, 3e7, 5e7, and le8 years. The two systems in the upper right above the tracks are both 
unresolved subarcseco nd binaries: the o t her fo ur that all fall on the tracks are considered 
to be single, although iDoppmann et al.l (l2008l ) suggests that Haro 6-lOS is a spectroscopic 
binary. 
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Fig. 10. 
spectral 
feature. 



5 10 15 20 25 30 35 
Ay [L-band ice absorption] (mag) 

— Extinctions derived from the JHK color-color diagram (open squares) and X-band 
fitting (filled circles) as a function of extinction derived from fitting the 3.1 /im ice 
The dash-dot line indicates equal extinctions on both axes. 



